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Abstract 

For the periodic sphaleron potential in the electroweak theory, we find the one¬ 
dimensional time-independent Schrodinger equation with the Chern-Simons number 
as the coordinate, construct the Bloch wave function and determine the correspond¬ 
ing conducting (pass) band structure. We show that the baryon-lepton number 
violating processes can take place without the exponential tunneling suppression (at 
zero temperature) at energies around and above the barrier height (sphaleron en¬ 
ergy) at 9.0 TeV. Phenomenologically, probable detection of such processes at LHC 
is discussed. 
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1 Introduction 

It is by now well-known that both the baryon number (B) and the lepton numbers (L) are not 
conserved in the standard electroweak theory because of the presence of anomalies [l||4| and 
instantons [5]. The non-Abelian nature of a Yang-Mills theory leads to a topologically non-trivial 
vacuum structure. An instanton solution (with index N = 1) in the 4-dimensional Euclidean 
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Yang-Mills field equations yields a tunneling amplitude from one vacuum |n) to another vacuum 
\n + 1) of order exp(— 2n/a), where a = g 2 /(Aii) is the gauge coupling; for small couplings, 
such tunneling is typically exponentially suppressed. Since there are multiple vacua labeled by 
the number n (n = — oo,... — 1, 0, +1,... + oo), separated by the same barrier, such a periodic 
effective potential implies that, in the SU(3 ) QCD case, the actual QCD vacuum should be 
described by a Bloch wave function labeled by 0, i.e., |0) = Jj) n exp(m0) \n) |6,7|. Because of 
the scaling property of the instanton solutions and the running of the coupling olqcd hi QCD, 
it is difficult to determine the barrier potential in QCD. On the other hand, the electroweak 
theory has a natural scale, namely the Higgs vacuum expectation value v, or equivalently, the 
W-boson mass mw = gw/2, where g is the SU( 2) gauge coupling. The existence of the closely 
related sphaleron in the electroweak theory [8] was first studied by Manton |9j and Klinkhamer 
and Manton [l0|. Although the sphaleron is not a topological soliton, it does have a Chern- 
Simons (CS) number (half-integer) and is important to the dynamics in the electroweak theory. 
The sphaleron energy E sp h measures the height of the potential barrier to the baryon- and 
lepton-number (B + L)- violating processes (which conserve the (B — L) number). 


The sphaleron has been extensively studied because it is likely to play a crucial role in 
the matter-anti-matter asymmetry in our universe due to the electroweak phase transition in 
the early universe 11 , a subject that has been extensively studied (see e.g., Ref 12 16 for 
reviews). So it is obviously important to know whether such processes can be observed and 
studied in the laboratory today. The possibility of sphaleron mediated baryon and lepton 
number violating processes at high energy colliders has been studied to some extent |1T 21 


When the energy reach is much lower than the sphaleron energy E sph = 9 TeV, these baryon 
number violating processes are exponentially suppressed, by a factor like exp(— ‘in/aw) ~ 10 -162 
where aw — 1/30. When the parton-parton energy approaches E sp h, such processes are not as 
suppressed. However, the result of the analyses done for collider physics so far is somewhat 
ambiguous. 


Here we like to estimate the (R + L)-violating rates in colliders using the periodicity property 
of the sphaleron potential as well as the presence of the kinetic term for the CS number n, which 
is necessary for a full first quantized treatment of the problem. To our knowledge, the discrete 
symmetry has not been emphasized in any of the sphaleron studies so far. As is well known, 
the quantum properties of a periodic potential is quite different from that of a single potential 
barrier. In this paper, we find the corresponding Bloch wave function and the conducting (pass) 
band (one-dimensional Brillouin zone) structure in the electroweak theory. 


Although n always takes an integer value at a vacuum state and a half-integer value at a 
peak of the potential (i.e., sphaleron solution), it takes continuous values as we move from one 
vacuum state over the sphaleron potential barrier to the next vacuum state. Since the sphaleron 
potential has been calculated already |b22|, we need only to evaluate the “mass” m in the kinetic 
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term rh{n)h 2 / 2 (dot denotes time derivative) to obtain the corresponding one-dimensional time- 
independent Schrodinger equation in which the CS number n is the coordinate to be quantized. 
For values away from half-integers, n 7 ^ 0, ±1/2, ±1,the definition of n is not unique, as 
the CS current is not gauge-invariant. We find that / 1 / 7 r, instead of the standard choice of n, 
is the most appropriate CS number that takes continuous values (where 7iir = /i — sin( 2 /i)/ 2 ). 
Introducing Q = /j/mw (so Q has the dimension of a coordinate) we obtain a constant mass m 
and 

{-Lw + v{Q ^ {Q) = mQl (L1) 

Using the known Higgs vacuum expectation value v = 246 GeV, W Boson mass rriw = 80 GeV 
and the Higgs Boson mass rri h = 125 GeV, we obtain 

V(Q) — 4.75 TeV (l.31 sin 2 {mwQ) ± 0.60 sin 4 (mvi/<5)) , 

E sph = max[V(Q)] = V = 9-11 TeV, 

m = 17.1 TeV, (1.2) 

where the potential V(Q) was obtained by Manton (see Fig. [I]). Determining the value of this 
mass m is a main result of this paper. Note that a rescaling of Q rescales m, though the physics 
is unchanged. 

E(TeV) 



Figure 1: The periodic sphaleron potential V(Q) as a function of the coordinate Q in the 
electroweak theory. The barrier height is 9 TeV. The dimensionless /i = mwQ is related to the 
Chern-Simons number n via n = /x/7r — sin(2//)/(27r). The extrema of V(Q) are at sin(2/i) = 0: 
the minima (vacua) are at integers n = p /n = ..., —2, —1, 0, ±1, ±2,... and the peaks (i.e., the 
sphaleron) are at n ± 1/2 |9|. 


To hnd the mass m, we have to make a couple of appropriate changes (gauge rotations) in 
the existing works. This is necessary because, although the static potential is gauge-invariant, 
different choices of static gauges tend to yield different masses in the time dependent kinetic 
term. The static sphaleron potential barrier has been calculated in two ways, namely the Manton 
method [9] and the method due to Akiba, Kikuchi and Yanagida (AKY) 1221. After the necessary 
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modifications just mentioned, we find that the mass m in the Manton approach is a constant, 
as given in Eq.(1.2), while the AKY mass m(ri) diverges (close to a linear divergence) as n —> 0. 
This divergent behavior is close to the simple example of ij 2 /(A\y\) + \y\ —» x 2 + x 2 if y = x 2 . 
After the redefinition, the constant AKY mass m = 22.5 TeV (with co-ordinate Q) is somewhat 
larger than the mass m = 17.1 TeV in the Manton case, while the potentials are close but not 
the same. So the overall features stay the same. Clearly a fully time-dependent evaluation of m 
may improve its value. Fortunately, the present approximation is good enough for our purpose. 

It was pointed out in (9,23,24 that turning on the U( 1) coupling (i.e., Weinberg angle 
sin 2 6 W = 0.23) will lower the sphaleron energy by about a percent. So it is reasonable to use 
E sp h = 9.0 TeV in phenomenological studies. 


Once we have the one-dimensional time-independent Schrodinger equation (1.1) with the 
mass m and the periodic potential V(Q) (1.2), it is straightforward to solve for the Bloch wave 
function, the conducting (pass) bands, their widths and the gaps between the bands. In Table 
[lj we give the lowest few bands and the ones that are close to the barrier height E sph = 9.11 
TeV. (Because of the higher mass m and higher potential away from the extrema, there are more 
bands in the AKY estimate.) We see that the first band occurs at about 35 GeV. With energies 
inside a pass band, the wave function spreads across the whole potential and transmission from 
one vacuum to another (at different integer n) is no longer tunneling suppressed. However, the 
width of the bands at low energies are exponentially narrow. Averaging over a few bands and 
their gaps at low energies, we find that the probability to lie inside a band is exponentially 
suppressed. This is simply another way to see the tunneling suppression effect. As the energy 
approaches E sph from below, the widths of the bands become bigger, while the gaps become 
smaller, so the effect of the pass bands becomes important. This is when the (B + L)-violating 
processes are no longer tunneling suppressed, even when the energy is still a little below the 
barrier height. 


It is most important to search for these (B + L)-violating processes in the laboratory. Par¬ 
ticularly interesting parton (left-handed quarks) scatterings in proton-proton collisions are the 
An = — 1 quark-quark annihilations at close to or above 9 TeV, 

qL + qL -> UfJrqqqqqqq + X, (1.3) 

which results in 3 anti-leptons (one from each family) and 7 anti-quarks plus other particles. 
The preferred quark content should contain 3 anti-quarks from the second family (c or s) and 
3 anti-quarks from the third family (7 or b). The above scattering is an inclusive process, so 
X (with net B = L = 0) may include any number of W ± , Z and Higgs bosons, mesons and 
photons as well as fermion-anti-fermion pairs. X has a net electric charge to maintain charge 
conservation of the process. These events probably look like fireballs. Similarly, one can consider 
a particular An = +1 quark-quark scattering, 


qL + qL —>■ e /i r bbbcccddduu + A" 


(1.4) 
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So a single (£>+L)-violating event can produce 3 positive sign leptons plus 3 6-quarks (An = — 1), 
or 3 negative sign leptons plus 3 6-quarks (An = +1). We expect roughly equal numbers of 
events for each. A very crude order of magnitude estimate gives 10 4±2 such events in the coming 
Large Hadron Collider (LHC) run at 14 TeV proton-proton collisions. Because the fraction of 
quark-quark scatterings with energies close to or above the sphaleron energy E sp h = 9 TeV is 
substantially bigger at the 14 TeV than at the 13 TeV, detection may be easier by comparing 
the two data samples. If the LHC energy can be raised by a few TeV, the (B + L)-violating 
event rate should increase by more than an order of magnitude. 

The rest of the paper is organized as follows. In Sec. 2, we review the properties of the 
sphaleron and the barrier potential first constructed by Manton |9|, where the necessary modifi¬ 
cations for our purpose are already incorporated. In Sec. 3, we construct the Lagrangian and the 
Hamiltonian for the Chern-Simons number after finding its kinetic term. After quantization, we 
obtain the one dimensional Schrodinger equation with Q = /i/m-w (instead of the Chern-Simons 
number) as the coordinate. It is then straightforward to find the pass band structure in the 
electroweak theory, which we present in Sec. 6. In Sec. 4, we point out the differences we have 
to make in the Manton construction to obtain what we believe to be the correct mass. In Sec. 5, 
we review the AKY construction of the potential and find the corresponding mass m(n), which 
diverges as n —> 0. After a transformation to the canonical form, we find that the AKY potential 
is not as different from the Manton potential as one might initially be led to believe. Sec. 4 and 
Sec. 5 may be skipped without loss of continuity. In Sec. 7, we discuss some phenomenology 
of the (B + L)-violating processes. Our main point is the possibility that they may be detected 
at LHC in the near future. Clearly more theoretical as well as phenomenological studies are 
warranted. Sec. 8 contains the summary and some remarks. Appendix A contains some details. 


2 Review 

Let us review the basic facts about anomalies, instantons, tunneling and the sphaleron potential. 
We make some changes to the original work in the study of the sphaleron as this will turn out 
to be important later. 

2.1 Background 

We start with the 577(2) weak interaction gauge fields A“(x) coupled to a doublet Higgs field 
$(x) and left-handed fermion doublets = (g{’“,/{), where / = 1,2,3 is the family index 
and a = 1, 2,3 is the color index for the quarks. To simplify the discussions, we ignore the 6/(1) 
gauge field (i.e., set Weinberg angle 6\y = 0) and the right-handed fermions; for our study, a 
non-zero 0\y as well as the Higgs coupling to right-handed fermions introducing fermion masses 
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and Yukawa couplings will have only minor changes to the analysis. We shall comment on these 
minor complications later. 

Consider only the simplified version of the standard electroweak model (with ds 2 = dt 2 — dx 2 
and H — c— 1), 

C = - 1 - Tr [F,j V F^\ + ^ ^ (d^h - v 2 ) 2 + (2.1) 

where 

(j a 

Fjw = F^— = d^A v d v An ig [A^, A u ], 

= dfl - igA 

- igA^f, (2.2) 

where A^(x) = A°cr a /2 and i = 1, 2,..., 12 for the 12 doublets of left-handed fermions. 

At the classical level, there exist til = 12 (i — 1, 2,..., til) globally conserved 17(1) currents 

jhY = 


corresponding to the conservation of the fermion numbers. However, this conservation is broken 
by the presence of anomaly jT, 21, 


r) /(*)/* — JL _ Tr 

M 16vr 2 




where F^ v is the dual of F flv . In the presence of instanton solutions in Euclidean space-time [5], 

_2 r 


N = 


9 


167T 2 


d A x Tr 


F^ v 


(2.3) 


where the topological index N takes only integer values. An instanton with value N leads to 
the tunneling process | n) —> \n + N). As the change of fermion number AAT 1 ' is the same for 
each doublet, we have the change in the lepton number L given by 


AN e = AN^ = A N r = N, (2.4) 

and the change in the baryon number B given by 

AB = ^ (3) (3)77 = 377, (2.5) 

o 

since each quark has B — 1/3 and there are 3 families and 3 colors |3|. As a result, the (B — L ) 
number is conserved, as AB — AL = 3 N — 3 N = 0. Since the electric charge Q e is always 
conserved, we have 

A(B + L) = 6N, A(B — L) — AQ e = 0. (2.6) 
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For example, an e + e~ collision can produce 3 baryons : e + e _ —y e~p~r~B+B+B+ where each 
baryon B+ is a bottom-charm baryon which decays to a nucleon plus mesons, e.g., B+(bcu) —> 
p7r + 7r _ 7r°. By comparison, a process in a proton-proton collision pp —> e + p + v T n can be mediated 
by an instanton also, though this process is further suppressed by the CKM mixing. At low 
energies, since an instanton action S = 27 t / aw where the weak coupling aw ~ 1/29.7, the 
tunneling rate goes like exp(— 2S) ~ 10 -162 , which is totally unobservable. 


In this paper, we propose that, under appropriate conditions, (B + L)-violating processes 
may not be exponentially suppressed. It is known that such (B + L)-violating processes become 
much less suppressed at high temperatures 25-27 and it may be responsible for the matter- 


anti-matter asymmetry of our universe [11]. Here we restrict ourselves to the zero temperature 
case and such processes in the laboratory only. 


The SU( 2) gauge theory has Euclidean instanton solutions as the result of mapping spa¬ 
tial 3-sphere S 3 (in 4-dimensional Euclidean space) to the gauge manifold, with homotopy 
7r 3 (SU( 2)) = 7t 3 (S' 3 ) = Z. That is, the Euclidean equations of motion (with t as the Euclidean 
time) have solitonic solutions where a spatial 3-sphere maps to the gauge held Ai(t,r,9,p). 
In Minkowski (3+1) space in the electroweak theory, one considers static Ai(p,r,9,p) and 
<h(/i,r, 9, p) where 0 < p < n. (We choose \n = 0) = \p = 0) as the reference vacuum.) At fixed 
r, the point p(p, 9 , ip) spans a 3-sphere, as shown in Fig. [2j The S 2 swept by the usual polar 
angles ( 9 , p) has zero size at p = rnr = 0 and p = mt = n, and maximum size at /i = nir = 7t/ 2 
corresponding to the sphaleron, an extremum (but unstable) solution of the static equations of 
motion. Time-dependence is then introduced when we treat p(t) as a function of time. 



Figure 2: Here the 3-sphere is spanned by the point p(/i, 9, p) where 0 < p < tt and the usual 
polar angles (9,p) spans a 2-sphere. The sphaleron corresponds to p = tt/2, when the 2-sphere 
attains maximum size. This unstable 2-sphere shrinks to a point at either vacuum (p — 0 or 7r). 
There are two ways to cover the 3-sphere : (a) the Manton (or the hnite mass) path (4.3), (b) 
the constant mass (or the spherically symmetric) path (2.7) used in this paper. 


A sphaleron is an unstable static solution of the classical equations of motion of the elec- 
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troweak theory (9,10 so its energy measures the height of the potential barrier that separates 
the | n) vacuum from the \n + 1) vacuum. Here we are interested in the potential between 
two adjacent vacua, not just the height of the barrier. There are two approaches in finding this 
sphaleron potential, namely the Manton approach [9] and the AKY approach 1221. Let us review 
these two approaches and comment on their relation. 


2.2 The Sphaleron Potential 

It is convenient to work in the temporal gauge A 0 = 0 and then further impose the polar gauge 
rA r = 0 to fix all local gauge freedom. This implies that A r = 0 for all r except may be at 
r = 0. In fact, as we shall see, A r = —/i5{r)f ■ a/2g. Here, we can ignore this for the moment. 
In the spherically symmetric ansatz, we can write the fields in the following forms 


$ = v (1 — h{r)) U I I + h{r 
cos [1 

Ai = -(1 - f(r))UdiU\ 

9 



U = 


r—»o r 


cos // + i sin fi cos 6 — sin /i sin 8e lv> 

sin fjL sin cos fi — i sin /i cos 6 y 

= h( 0) = 0, /(oo) = h( oo) = 1, 


(2.7) 


where {/j,,9,<p} are the polar angles of S 3 and the SU( 2) matrix U maps this S 3 to the group 
manifold Sg Lr ^ 2 y With the above Higgs held profile h{r ) and gauge held prohle /(r), /r = 0 
and /d = | correspond to the vacuum and the sphaleron respectively. Note that the asymptotic 
Higgs vacuum expectation value is always = [0, u] for every 0 < fi < n (shown in Fig. [3] (b)), 
which is different from the Manton setup (9) (shown in Fig. [3](a) ). This difference is explained 
in Sec. [4] as it is important in the determination of the mass. To avoid confusion, we shall call 


the above ansatz (2.7) the constant mass construction. 


The energy of the sphaleron path is given by 


47T 

Vm{9) = 


dr \ + wt/(/ “ !)] shr V + 2r h ' + 4m w(/ - h ) 

1 


+4 m^[f(h — 1 ){hf + / — 2 h)\ sin 2 n + -m H r (h — l) 2 sin 2 /i f sin^ /j, , (2.8) 

where the static equations of motion at the sphaleron (/i = 7 t/ 2) are 


r 2 f" = 2/(1 - /)(1 - 2/) + m 2 w r 2 h(f - 1), 
( r 2 h'y = 2h(l — f ) 2 + ^ m 2 H r 2 (h 2 — 1 )h. 


(2.9) 






Using the above boundary conditions, f(r) and h(r) can be easily solved numerically (see Ap¬ 
pendix A). As mentioned in Ref 10 , the sphaleron solution indeed has topological number Its 


energy E sp h = 9.11 TeV measures the potential barrier height. As shown in Fig. [2](b) , varying fi 
from 0 to Ti spans a 3-sphere, which also goes from the |n = 0) vacuum over the potential barrier 
to the \n = 1) vauum. So it is a reasonable approximation to take the /(r) and h(r) solution of 
Eq.(2.9) for /i = 7t/2 and insert them directly into Vm{p) (2-8) to obtain V(fx) (1.2). 


3 The Effective Time-independent Schrodinger Equation 
for the Chern-Simons Number 


In this section, we write down the effective one-dimensional time-independent Schrodinger equa¬ 
tion where the coordinate is the CS number n. Since the potential V(n) (1.2) has already been 
evaluated, we only need to find the mass m. 


The Lagrangian (2.1) has kinetic terms for both A p {x) and $(a;) that include the time 
derivatives. Here, we are only interested in the change of the CS number n (or p) as a function 
of time : we start with the static solutions for Aj(?r, r, 6 , (p) and <h(?7, r, 6, ip) and then introduce 


time-dependence only from the change in n, i.e., n — 


dAi dAi dn 


dt 


dn dt ’ 


< 9 $ 

~dt 


<9<L dn 
dn dt' 


So the kinetic terms in the Lagrangian (2.1) yield a kinetic term for n. That is, we set to 
determine the mass (which can be a function of n) to obtain the kinetic term m(n)h 2 /2. Once 
we have the Lagrangian for n, we quantize the system and write down the time-independent 
Schrodinger Equation for n. 


A comment is in order here. In principle, we should allow explicit time dependence in both 
A p (x) and <h(x) to obtain the kinetic term for /i in a fully gauge invariant approach. This 
requires allowing A 0 to be non-zero. A 0 then is determined implicitly by imposing Gauss law, 
the held equation obtained by varying with respect to A 0 (or perhaps A p ). We can still impose 
other gauge conditions, for example, spherical symmetry together with a radial gauge condition 
rA r = 0, so the mass density should be spherically symmetric. In practice, the role of Aq is 
to subtract off from held variations the part that is just an inhnitesimal gauge transformation. 
As a result, our estimate ignoring A 0 (or setting it to zero) should give an overestimate of the 
kinetic energy and hence the mass parameter m. Fortunately, as we shall see, this over-estimate 
will have little or no effect in the phenomenology discussed in this paper. 


It turns out that the variable p used in Sec. 2 is a more suitable coordinate for our purpose. 
It is closely related to n (this relation will be explained in Sec. 4), 

sin(2/i) 


rnr = /i — 


(3,1) 
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which equals /i = nn at integer and half-integer values of n. i.e., solutions (extrema) of the 
static equations of motion. Let us promote /i to a time-dependent variable fi(t) and compute 
the Lagrangian in terms of {/q /i} 


L= d 3 x£ = -m 


fr 


2 {m w y 




where the constant mass m is given by 
8TT7TIW 


m = 


r 



r) — l) 2 + 2 (m w r) 2 (h(r) — l) 2 ] = 17.1 TeV., 


(3.2) 


(3.3) 


where f(r) and h(r) are the /j = 7t/2 solutions obtained from Eq.(2.9). Define a coordinate 
with length dimension Q = /i/mw such that L = \mQ 2 — V(Q) (to make the commutator 
dimensionless as well). From this Lagrangian, one can define the canonical conjugate momentum 
of Q, and hence the Hamiltonian H , 

(3.4) 


dL 


vr, 


7r( 3 = S = m< ^’ H = nQ - L = ^+ V (Q), 
dQ 2 m 

By imposing quantization on the variable Q, [Q,tiq\ = i, ttq in H can be replaced by —^ i-^, 
resulting in the following familiar one-dimensional time-independent Schrodinger equation, 

1 3 2 


2 md Q’ + V ^) »«) = «*«). 


(3.5) 


where V L (Q) is the eigen-wavefunction of Q and energy E is the eigenvalue. This equation is 
solved in Sec. [6} Readers mostly interested in phenomenology may go directly there. 


4 Comparison to the Original Manton Construction 

As pointed out earlier, the above evaluation of the fi -independent mass m is obtained after 
making two changes to the original Manton construction. Here we like to explain these two 
changes. Both are gauge changes so the potential V(fj) remains the same as that obtained by 
Manton. However, the static gauge transformations are no longer pure gauge changes when we 
attempt to extract the time-dependent kinetic term from the static solutions. The differences 
are illustrated in Fig. [2] and Fig. [3] and summarized in the end of this section. 

4.1 Rendering the Mass Finite 

The original Manton construction starts with (following his notations) the following ansatz 
involving the unitary SU( 2) matrix U°°, 
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®(v,r,0,<p) = v[l~h{r)\\ _.° ) + vh(r)U°°{n, 9, cp) 1 ° 


e cos /i 


1 ’ 


^ i (^r,0,^) = -/(r)i/ oo 9 i (Z7 oo t ) 

9 


U°° = - 


1 / 4>f 


e J/i (cos )jl — i sin /i cos 6 ) 


sin fi sin 9e l ' fi 


$°°(jM,¥>) = v 

f(r) 


v J \ — sin/isin 6e e * M (cos/r + isincos#) 

sin fi sin 9e lv 


e lfl (cos fi + i sin fi cos 9) 


lim —— = hi 0) = 0, /(oo) = hioo) = 1. 

r— >o r 


(4.1) 


As fi ranges from 0 —* it, the Higgs field at r —» oo takes a hedgehog form, as illustrated in 
Fig. | a. Note that this choice is not suitable for onr purpose since the asymptotic Higgs field 
is in different gauges as fi goes from 0 to §. It results in a divergent mass m. This is because 
a change in /j results in a change in $ at r —» oo, even when the sphaleron is supposed to be 


localized. To remove this divergence, we simply make a gauge change to the fields in Eq. (4.1) 
so at r —> oo, —» [0, v] independent of /i. That is, 


$ = = v (i _ ^[/oof [ 0 ] + v h 1 ° 

\e lfl cos fi J 

Ai = U^AiU 00 + -U^diU 00 = -(1 - f (r))U°°^diU° 

9 9 

^ I e _l/i (cos /j, + i sin fi cos 9) — sin /i sin 9e lv 

l sin /i sin 9e~ l>p e Zfl (cos /i — i sin /i cos 9) 

lim^—- = hi 0) = 0, /(oo) = hioo) = 1. 

r->0 r 

This renders the mass m finite so we refer to this as the finite mass construction. 


(4.2) 


4.2 Maintaining Spherical Symmetry 


Note the e ±l,J ‘ factor difference between Eq.(4.2) and Eq.(2.7). Let us write the unitary matrix 


f/°°t in Eq.(4.2) in terms of the Pauli matrices and a 4-component vector x, 

U °°* = x 4 I + i(xjCTi) 
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Figure 3: The direction of the vacuum expectation value of $ at large distances. [Re $!, Re <3> 2 ] T 
in (a) hedgehog gauge ( |4.1 ), where <f> is rotated as /i varies at asymptotic distances, to reach a 
hedgehog shape at /i = 7 t/ 2; this introduces an infinite kinetic term; (b) unitary gauge at the 
sphaleron (2.7), where <f> always stays the same at large distances even as /i varies. 


p(v,o,<p) = 


(4.3) 


so that its determinant yields x • x = — 1- We associate this unit vector x with a point 

p(/i, 9, ip) in S 3 , 

( — sin /i sin 9 sin ip \ 

— sin /i sin 9 cos tp 
sin fjL cos /i(l — cos 9) 

\ sin 2 /i cos 9 + cos 2 /j / 

This spans the S 3 as shown in Fig. |2](a). On the other hand, U in Eq.(2.7) yields 

/— sin /j sin 9 sin ip\ 

— sin /i sin 9 cos <p 
sin fi cos 6 

\ COS /l J 

which spans the same S 3 but in a different way, as shown in Fig. |2](b). 

As a point p(/q 9, ip) in S 3 , each of the above two cases generates the same S 3 by generating 
S 2 for fixed // and then varying 0 < p < 7r, but via two different paths. The mass m from the 
finite mass path (4.2) yields 


p{p,o,<p) = 


(4.4) 


m(pi) 


47T m w 

g 2 


87im w 

g 2 


d(cos9)d(mw'r)l [4 + 2(1 — 4cos0 + cos 2 9) sin 2 /i] (/ — l) 2 
+ ^( m w r ) 2 [4 + (1 - cos 2 9) sin 2 2/x] (/i — 1) 2 |, 
d{m w r) [4(/ - 1) 2 F(/i) + 2(m^r) 2 (/i - 1 ) 2 G(/i)] , 


(4.5) 
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where 


(4.6) 


2 1 

F( n) = l+o sin2 Ah C(/-0 = 1 + « sin2 2 /b 
3 0 


while our path (2.7) yields m in Eq.(3.3). We see in the mass integral (4.5) in the finite mass 


path (4.2) that the ^-dependent terms in m(n) are introduced by the angular dependent terms, 


which breaks the spherical symmetry adopted in the original ansatz. Observe in Fig. |2](a) that, 
as /i varies, p(p, 0 = 0, ip) is fixed at (0, 0, 0, l) r while p(/i, 0 ^ 0) travels at different “velocity” at 
different 6. This 0-dependent motion in S 3 introduces the 0-dependence in the mass m density, 
leading to a bigger mass m{fi) that depends on /i. To maintain the spherical symmetry, we 


instead choose U in (2.7, 4.4); that is, we simply drop the e ±l,i factors in Eq.(4.2). This yields 


the constant mass m ( 3.3| ) (i.e., with F(n) = G(/i) = 1 instead of the values given in Eq.(4.5)). 
Note that the difference between these two cases can be treated as a pure gauge transformation, 
so the static potential H(p) remains the same in each case. However, it is no longer a simple 
static gauge choice when we attempt to extract the time-dependent kinetic term from the static 
solution in the static gauge. 


4.3 Relation between the Chern-Simons Number and the Canonical 
Variable p 

To discuss the topological baryon number, the Chern-Simons number is often used. However, 
one has to pay extra attention to it as it is not fully gauge invariant in a time slice and it is not 


necessarily equivalent to the topological baryon number. Writing the Chern character (2.3) in 
the Chern-Simons form, 


N = 


9 


167T 2 


iW = 


d A x TV 




32vr 2 


F^ 


F a A a 

Vp-L x O 


= / d^xdnK^, 


9 abc a a Ab a c 


So 


N = / d 3 xK° I t + 

/ \t=to 


rto 


K ■ dS. 


(4.7) 


(4,8) 


— CXD J S 


We choose ji(t = — cx)) = 0 and /i(t = t 0 ) = p. 


Let us start with the fields in Eq.(4.1), but dropping all the exp(±f/i) factors, as we have 


just discussed. Since this (as well as the A t in Eq.(4.1)) drops off like 1/r asymptotically, 


the surface term in Eq.(4.8) does not vanish. One can make use of the residual 77(1) gauge 

-iQ(r)f ■ a ), 72(r) = /ztanh (/3r), 


transformation 28 


exp 


(where j3 is large) to rotate A t to A t in Eq.(2.7) at large distances. Now, A t drops off exponen¬ 


tially so the surface term no longer contributes to N. This difference is related to the asymptotic 
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direction of the Higgs vacuum expectation value (see Fig. [3]). Such a transformation contributes 


a term to the Chern-Simons number [29], so the 4-dimensional Chern number (2.3) becomes 

2Q(r) — sin(2Q(r)) r =°° 2/j, — sin(2//) 


N = 


d xK \t=t 0 + 


2vr 


r=0 


2vr 


(4.9) 


since K° vanishes. This is the relation (3.1). To get to the Ai at all r as given in Eq.(2.7), we 


have to take /3 —> oo. In this case, we find that A r = —fi5(r)f ■ <3/2g instead of A r = 0. (Note 
that the polar gauge condition rA r = 0 is still satisfied.) This 8(r) term has no impact in the 
evaluation of the potential V (//) or the mass m. 

Let us summarize in the diagram below the changes carried out in this paper. 


Manton (4.1) 


Fig. 3(a) ->• 3(b) 


Finite mass (4.2, 4.5) 


(mass density violat¬ 
ing spherically sym¬ 
metry ) 


Fig. 2(a) 2(b) 


remove e 


=t ifi 


Fig. 2(a) 2(b) 


f K ■ dS ^ 0 


AKY 1221 


Fig. 3(a) —>■ 3(b) 


Constant mass (2.7) 


Sec. [5] 


(spherically symmet¬ 
ric mass density) 


5 Comparison of the Manton and the AKY Construc¬ 
tions 

Here we compare the sphaleron potential obtained in the Manton approach versus that in the 
AKY approach. This comparison is possible only when both use a coordinate with a canonical 
kinetic term. We see that their masses differ by about 30 % while the potentials are quite close. 

5.1 The AKY Potential for the Sphaleron 

The Manton construction of the path between a sphaleron and the vacuum comes from the 
picture of shrinking a S 2 loop in S 3 to a point. The study of this path is performed in a 
different manner by Akiba, Kikuchi and Yanagida (AKY) 1221, in which a static minimum- 
energy path is obtained by considering the most general spherically symmetric ansatz of the 
fields, 
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Aq =-a 0 (r,t)x a , Aj — - 
9 9 


ai(r,t)xjx a + 


f A (r,t) - 1 


£jam%m H~ 


/bM) 


(Sj a XjXa) 


$ = (. H(r , t ) + K(r, t)ia ■ x) 


: 


(5.1) 


It is shown in 28,30 that under these spherically symmetric held configurations, the SU( 2) 
gauge and Higgs sectors in the model (2.1) reduce to a (1 + l)-dimensional model, 

1 ..„ 1 


& = ~ -/-U„r + (D.xYD-x - - l ) 2 


2m 


w 


r 2 (D^Dy - -(|x| 2 + 1)|0| 2 + Re (x*0 2 


m 


H 2 


r‘m‘ - 1)' 


(5.2) 


where, 


ds 2 = dt 2 - dr 2 , f ltu = d^a v - c0a /t , x = /a + h/s, <j> = H + iK , 

'l 

D„x = ( 3 „ - ia„)x , = (d p - -a p )</>. 


(5.3) 


It preserves the 17(1) subgroup of 577(2) via the gauge transformation + C0I2, x ~^ 

exp(iO)x and 0 —> exp (|0) 0. In a time-independent setup, a 0 = 0 and cp = 0 are chosen to 
fix the local gauge. However, static equations of motion from (5.2) have no non-trivial solution 
when the winding number N is not a half-integer, since they are not extrema of the energy 
functional. The AKY static path is solved by minimizing the energy functional with respect to 
a fixed N (2.3), namely, minimizing the functional, 

W\x, 0 ] = E stat [x , 0 ] + v(N[x, 0 ] - n), 

where r] is a Lagrange multiplier. The static equations of motion become 


(5.4) 


-dr + ~M - 1 )+m 2 w\ 


X = m w <j> - imwCdrX, 


~dr + ^(l^l 2 + 1) + “ 1) 


( r 0) = -x0* 


(5.5) 


with ( = awv/2'nmw The asymptotic behavior of the solution is, 
r —>■ 0, 

X « e~ iq [1 + c A (m w r) 2 + i\{-c H c K + c A ()(m w r) 3 ] 
0 « e _ *2 [ Ch + -Ac h (c 2 h - l)(m^r) 2 + ic k (m w r)\ 

r —> oo, 

X ~ 1 + Re ( d A e~ ar ) + i Re (d B e _Qr ) 

*Re {^e~ ar ) ’ 


2 -j_ dH c ~»rinr 


(5.6) 


m w r 
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where 


a = 


m w 


(y4-c 2 ±;c), d 


B — 


C rriwOtdA 

™w - a 2 ' 


(5.7) 


Note that a different boundary condition is chosen here compared to that in Ref 22 . These 


two sets of boundary conditions are equivalent up to a rigid U{ 1) transformation of hi = —q, 
where q is constant everywhere. Likewise, due to the aforementioned reason (see the summary 
diagram at the end of Sec. [4]), one has to keep the large distance vacuum in the same gauge in 
order to discuss the kinetic term. 


Since the original AKY set of boundary conditions chosen in 22 is not in the unitary gauge, 


the surface term f K ■ dS does not vanish; one has to make use of the residual U(l) gauge 
transformation with Q(r —>■ oo) —>• — q exponentially and fi(r = 0) = 0 to make A“ drop faster 
than 4 such that f K ■ dS = 0 |9|. Such a transformation contributes a term (q — sing)/27r to 


the CS number [291. The 4-dimensional Chern number (2.3) thus reduces to 


N = / d 3 xK° + 


q — sin q 
27T 


1 

27T 


dr Re (■ ix*d r X ) + Sb |^° + 


q — sm q 
27T ’ 


(5.8) 


where a reference vacuum of CS number 0 is chosen at past infinity. We make a gauge transfor¬ 
mation of fl(r) = — q lim^oo tanh(/3r) to the AKY boundary conditions to obtain our boundary 


conditions in order to preserve regularity at the origin. This changes ai(r) in (5.1) to a 5- 
function. This only changes the CS number and its effect is already explicitly included in the 
term (q — smq)/2n while J d 3 xK° is safely calculated from the rest. 


The equations of motion (5.5) are solved numerically under the above boundary conditions 

fa). 


and the AKY potential (22 is reproduced, as shown in Fig. [4|1 


5.2 Approximate Equivalence of the modified and the AKY Con¬ 
structions 

Following the discussion in Sec. [3j we also examine the Lagrangian for AKY case. Promoting n 
to a time-dependent function n(t), the Lagrangian is written in {n,n}, L = | M(n)n 2 — V(n), 
where [] 


47T f 

M(n = — dr 

9 1 J 


(' dnf A ) 2 + ( d n f B ) 2 + 2 m 2 w r 2 (( d n H ) 2 + ( d n K ) 2 ) 


(5.9) 


Following the definition (5.4), one can verify by a Legendre transform that the potential 


V{n) has a linear cusp near integer n (see |22|), 


dE(n) 


dn 


n Anniw 


n=0 \n\ a w 


(5.10) 


l M(ri) has dimension [mass] 1 since n is dimensionless. 
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-1 


-1 


V(n) 



M{n) 




(a) 


(*) 



Figure 4: (a) Solid line: the AKY potential V{n ) is periodic in n which has a cusp shape 
potential around each topological vacuum. Dashed curve: M{n ) is divergent at integer n = 0,1. 
(b) M[n ) in log-log scale. It is approaching a straight line which gives M{n ) oc nr 098 as n —» 0. 


It is interesting to see from Fig. [4] that M{n ) diverges at n = 0,1. 

This divergence is actually a nice feature which helps to remove the cusp of the AKY potential 
around a vacuum. The simple example illustrates how this happens : y 2 /(A\y\) + \y\ —> x 2 + x 2 
by the reparametrization y = x 2 . Both the divergence in the mass l/(4|y|) and the cusp feature 
of the potential \y\ at y = 0 vanish in the canonical variable x. If M{n ) diverges as n -1 , AKY 
potential will become quadratic-like around vacuum which is similar to the Manton potential. 
We checked numerically the power of divergence in M{n). It approaches nr 098 at n = 0.01, 
instead of n -1 , though the latter value is possible as n —>■ 0. For phenomenological purposes, the 
exact limit of divergence power at very small n does not affect our main result. We extrapolate 
the divergence power as n~°' 9S for small n in discussions below. A similar behavior has been 


observed in an instanton calculation 31 


After a reparametrization y, = a f \J M[n)dn , where a is fixed by setting // = | at the 
sphaleron, the Lagrangian with a canonical kinetic term is obtained, L = m/i 2 /(2m^/) — V(y). 
Numerically, m ~ 22.5 TeV which is around 1.3 times the mass m = 17.1 TeV in (1.2). We 
observe that the shape of AKY potential V(y) displayed in Fig. [5] is now quite close to the 
Manton potential in this variable. 


We again quantize the variable Q = y/mw and solve for the eigenvalues of the Hamiltonian 
to get the band structure for the AKY case. Results are presented in Sec. [6j Based on our 
analysis, we see that y/n (instead of n ) is the best choice as the CS number. 
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Figure 5: Comparison of AKY potential (solid line) and Manton potential (dashed line), (a) in 
CS number n, (b) in cannonical variable //. In the Manton case, n = (2ji — sin2/x)/(27r), which 
reduces to n = /x/7r at half-interger values of n. 


6 Bloch Wave and Band Structure 


6.1 Solving the Schrodinger Equation 


In this section we solve the Schrodinger equation (3.5) in both the Manton and the AKY 
construction. Since V ( Q ) is periodic, the energy eigen-function of Htjj = Ei/j should be Bloch- 
waves, 2 

( 6 . 1 ) 


(<21*.') = i>k(Q) = e’* Q m(<2), u k (Q) = u t I <2 + 


7T 


m w 


It is well known that the energy spectrum of a periodic potential has a band structure: continuous 
bands (solutions to the Schrodinger equation (3.5)) of certain widths separated by bandgaps (i.e., 
regions with no solution). 


It is standard to call the edge energies of the bands as eigenvalues and the corresponding 
wave functions as eigenfunctions. Note that V(Q) is symmetric with respect to both Q = 0, i.e., 
V(Q) = V(—Q), and Q = n/2mw, i.e., V(Q + w/2m w ) = V(—Q + n/2m w )] so there are 4 types 
of eigenfunctions, namely (SS), (S54), (AS 1 ) and (AA), where the first letter denotes symmetric 
(S) or anti-symmetric (A) about Q = 0 and the second letter denotes that about Q = / n/2mw- 
(SS) and (AA) eigenfunctions have period it/ mw while (SA) and (AS) eigenfunctions have 
period 2n/mw Starting from the (SS) ground state (i.e., lower edge energy of the first allowed 
band) as the energy is increased, the eigenfunctions have period [1,2, 2,1,1,2, 2,1,1,. (tt/vtiw 

2 Similar to the 0-vacuum in QCD, the Hamiltonian eigenstates should be superpositions of all local states, 
\k) = J2 n e ink ”/ mw | n). 
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E(TeV) 


E(TeV) 



9.20 


9.15 


9.10 


k(m w ) 


9.05 


(b) 


Figure 6: The (reduced) Brillouiu zone plot (6 conducting bands shown) near the barrier top for 


the Manton case: (a) The energy eigenvalue E(k) of the Schrodinger equation (3.5) is shown as 


a function of the propagation wave vector k (in units of m-w)', (b) The dark regions are the con¬ 
ducting bands while the gaps are regions where there is no solution to the Schrodinger equation 


((3.5) or (6.3)). The energy spectrum is almost continuous above the barrier while bandwidths 
below the barrier decrease exponentially as the energy goes down. The band structures in the 
AKY case are similar (see Table [T]) . 


One can approximate the band edge energies by the semi-classical method 32 , 
rQo 


.7T 


/ p(Q, E)dQ — K— ± arctan 
/ o 2 


tanh | / p(Q, E)dQ 
Q o 


, K e N 


( 6 , 2 ) 


where p(Q, E) = y/2m\E — V(Q)\ and Qo denotes the classical turning point in the interval 
[0, 7 t/2 mw\- The lowest band has edge energies with K = 0 ( SS ) and K = 1 (AS), the next 
band has edge energies with K = 1 (SA) and K = 2 (AA), and so on. This approximation is 


very accurate when the integral inside tanh in (6.2) is large 


ft is convenient to solve the Schrodinger equation (3.5) in momentum space using Bloch- 
waves, 


E 


-- (21 + k) 2 Sl m + Vl-r 

2m 


'U'k.l E Ufa r 


(6.3) 


where V)_ m , Uk, m are the Fourier coefficients of V(Q), Uk(Q) respectively. Here we solve Eq.(6.3) 
numerically for the band structure. Table [T| lists the band center energies and their bandwidths 
near the top of the potential barrier and around the bottom of the potential V(Q). Six bands 
near the top are shown in Fig. [6} In Fig. [7| we show that the logarithms of the bandwidths 
follows a linear curve. Since the band gap sizes change relatively slowly, it is easy to obtain the 
approximate band energies and their bandwidths for those not shown in Table [lj 

In the semi-classical approximation, a typical wavefunction with energy E will have an energy 
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Manton 

AKY 

Band Center 

Width 

Band Center 

Width 

Energy (TeV) 

(TeV) 

Energy(TeV) 

(TeV) 

9.113 

0.01555 

9.110 

0.01134 

9.081 

7.192X10" 3 

9.084 

4.957xl0 -3 

9.047 

2.621x 10 -3 

9.056 

1.718X10" 3 

9.010 

8.255x 10 -4 

9.026 

5.186X10" 4 

8.971 

2.382x 10~ 4 

8.994 

1.438xl0” 4 

8.931 

6.460x 10~ 5 

8.961 

3.747xl0 -5 

8.890 

1.666x 10~ 5 

8.927 

9.279xl0“ 6 

8.847 

4.114x 10“ 6 

8.892 

2.198X10" 6 

8.804 

9.779xl0~ 7 

8.857 

5.008xl0“ 7 

8.759 

2.245x 10~ 7 

8.802 

l.lOlxlO” 7 

8.714 

4.993X10" 8 

8.783 

2.341 xlO” 8 

8.668 

1.078xl0 -8 

8.745 

4.828xl0~ 9 

8.621 

2.262x 10 -9 

8.707 

9.673xl0 _1 ° 

8.574 

4.622x 10~ 10 

8.668 

1.886xl0 _1 ° 

8.526 

9.210x 10~ n 

8.628 

3.580xl0 -11 

8.477 

1.792xlO -11 

8.588 

6.622xl0 -12 

8.428 

3.411x 10~ 12 

8.548 

1.211xl0 -12 

8.379 

6.395x 10~ 13 

8.506 

2.167xl0 -13 

8.328 

1.208x 10~ 13 

8.465 

3.553xl0 -14 

0.3084 

- io- 169 

0.3146 

- IO” 204 

0.2398 

- 10- 171 

0.2454 

- 10~ 207 

0.1712 

- 10- 174 

0.1759 

~ io- 209 

0.1027 

- 10- 177 

0.1061 

- 10- 212 

0.03421 

~ 10- 18 ° 

0.03574 

- 10- 216 


Table 1: Some of the top and the bottom band energies and their widths (in TeVs) are shown. 
There are 148 bands up to E sph = 9.11 TeV in the Manton case and 164 bands in the AKY case. 
The band gap is about 70 GeV at low energies and decreases to about 30 GeV close to E sp h- 
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Width(TeV) 



E(TeV) 


Figure 7: Bandwidths of Manton path (solid line) and AKY path (dashed line) decrease ex¬ 
ponentially as energy decreases. This quantitative behavior obtained numerically fits well with 
that obtained in the WKB approximation (6.2). 


spread SE that is bigger than the width of band that it overlaps with. So it will have a tunneling 
probability amplitude that goes like 


A(E) ~ 


bandwidth 

bandgap + bandwidth 


where the bandgap is around 70 GeV towards the bottom of the potential. With an exponentially 
small bandwidth for low energies, the tunneling probability is clearly exponentially suppressed, 
as expected. As the energy E approaches E sph from below, the band widths increase quickly 
so the exponentially suppressing tunneling factor begins to vanish, so A (E) —> 1. As we go 
above E sp h, the band gaps shrink as well so essentially there is no tunneling suppression. This 
expected property is shown in Fig. [6] and Fig. [8j In fact, for any given k away from the edges, 
the wavefunction ipk(Q) for E(k ) > E sph is well approximated by a single free plane wave mode, 
i.e., we have the time-dependent wavefunction, 


k(Qj t ) ~ exp (ipQ - i£t ) 

where the ’’momentum” p and the effective energy £ depend on E{k ) and the average of the 
potential < V >= 4.197 TeV, 

2 

£ = E(k)— < V >= — 
v ; 2 m 

Note that the mean velocity is given by < v >= p/m, where \p\ > 13 TeV for E{k ) > E sp h. 
Note that the bandwidths A for bands close to the bottom are well approximated by tight- 
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binding model and the WKB method, 


A 


(n ± l\H\n) ~ exp 



(6.4) 


where c is a 0( 1) factor. Note that c = 1 for pure gauge theory, and c ~ 2.2 and c ~ 2.7 
respectively in the above two evaluations in Table [lj This can be understood as the size of the 
sphaleron is fixed by the electroweak breaking scale. Following Ref |4], we have 


A ~ exp[— Sa 


S h ] = exp 


2vr 

Ol]y 


—{m w p) 2 , 
O-w 


(6.5) 


where p is the sphaleron size, which is of order l/mw- This gives a value not far from the tight- 
binding result, where the band width measures the tunneling amplitude. It is not surprising 
that AKY has a bigger value for c, as its sphaleron size is in general bigger than that in the 
Manton case. In pure gauge theory, the instanton size p is integrated over so the second term 
in the exponent in Eq.(6.5) drops out. 


At very low two-body scattering energy E 2 (with some energy spread), we see that the 
probability of any (B + L )-violating process is exponentially suppressed. As we increase E 2 , 
such processes will be exponentially less suppressed. This behavior is shown by the solid curve 
in Fig. [8j The two pass band structures in Table 1 yield almost the same curve (which is really 
two curves overlapping with each other). As the energy passes E sp h, the exponential suppression 
factor vanishes. 


6.2 Comparing with Earlier Studies 


Earlier works are also presented in Fig. [8] for comparison. As we increase E 2) it is generally 
believed that such processes will be exponentially less suppressed. This behavior has shown up 
in perturbative calculations, and may be summarized in Fig. [8j where, for ease of comparison, 
we normalize all the curves for the various rates to the same value at E 2 ~ 0. They have been 


captured in the following formula 33-38 
<r(A n = ±1) ~ exp ( c 


47T 


Olw 


i + -(A 

8 \E„ 


4/3 


9 

16 


E 2 

E 0 


( 6 . 6 ) 


where E 0 = \/Qirmw / oiw = 15 TeV and c ~ 2. The perturbative analysis of | An| = 1 processes 
was started in [33i, 34], which show that the inclusive cross-section rises exponentially with 
energy. The authors pointed out that unitarity breaks down if one extrapolates the result to 


high energies. The E^ 3 term in the exponent factor (in square brackets) in Eq.(6.6) was later 

which includes only the E% 3 


obtained 


35 


39 


; this is captured by the dot-dash curve in Fig. 


term (and the —1 term) in the exponent factor (in square brackets) in Eq.(6.6). This clearly 
violates unitarity, implying that its validity is at best limited. Including the next El term in 


22 




















E(TeV) 


Figure 8: Comparison between the exponential suppression factor of (B+L)-violating processes 
in our case (solid line) and previous estimates. For ease of visualization, we normalize all cases 
to —1 at zero energy. The curves from the two band structures in Table 1 actually overlap 
to form a single solid curve that takes zero value above 9 TeV. The exponent factor (inside 
the square bracket) in Eq.(6.6) is plotted here. The first order (dot-dashed) and second order 
(dashed) approximation in Eq.(6.6) diverges differently at large E 2 . 


the exponent in Eq.(6.6), its behavior is shown by the dashed curve in Fig. [8j Additional 
terms with higher power dependence on E 2 are known to exist. So this approach yields an 
inconclusive result. Non-perturbative studies have also been carried out in toy models (see e.g., 
Ref [401). Such analyses seem to imply that the rate of (B + L)-violating processes is always 
exponentially suppressed to a level at which they will never be observed in the laboratory at 


any finite energies 13 


Expressed in terms of the WKB formula for tunneling through a potential barrier V(Q ) 
(with maximum value V max ~ mw/ctw ) at energy E ~ 0, one has a tunneling amplitude 
~ exp (— f \/2mVdQ). Comparing this to the tunneling amplitude ~ exp(— 2^/aw) estimated 
in earlier works in the semi-classical approximation, we find that mass m ~ m\y/&w- That 
is, the existence of the mass m and its order of magnitude value have been implicitly assumed 
already in earlier studies. We simply bring this semi-classical study to a fully first quantized 
analysis. Our approach includes the kinetic term for the CS number n/ tt, with Q = n/mw 
playing the role of the spatial coordinate. This allows us to carry out a fully first quantized 
analysis that incorporates the periodicity property of the sphaleron potential. The quantized 
version exhibits the important resonant tunneling effect via the Bloch theorem, with a simple 
band structure, where tunneling inside a pass band is unsuppressed. This crucial feature is 
absent in the earlier studies. 
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7 Phenomenology 


All the (B + L)-violating processes conserve electric charge, QCD color and the (B — L) number. 
The description of (B + L)-violating processes is cleanest in e + e - collisions, so we start with 
this. However, the energy needed will not be reached in the foreseeable future. So we turn 
to proton-proton collisions, which offer a much better chance of detection. This is particularly 
interesting at the Large Hadron Collider (LHC). 

7.1 e + e~ Collisions 

Let us first consider the following A(B + L) = 6 process in an e + e“ annihilation starting in the 
vacuum state |n = 0) going to |n = +1), 

e + + e —> l e + Ip + It + 9(i) + 9(i) + 9(i) + 9(2) + 9(2) + 9(2) + 9(3) + 9(3) + 9(3) + -V, (7.1) 

where the subscripts are family indices. Here, 12 fermions are produced, a lepton and 3 quarks 
from each family. This and all scatterings discussed in this paper are meant to be inclusive 
processes; that is, they may include (as denoted by X ) any number of W ± , Z and Higgs bosons, 
mesons and photons as well as fermion-anti-fermion pairs, where X has zero net baryon number 
and zero net lepton numbers but with an electric charge to preserve charge conservation of the 
process. With 3 b-quarks from the 3rd family, the 12 fermions require about 20 GeV of energy, 
so we do not expect any phase space constraint in the (B + L)-violating processes. We do not 
address the issue of the number distribution of the massive bosons in this study. Instead of 
3 quarks from each family, all quarks from the same family can also be produced, but such 
processes will be CKM mixing angle suppressed. Equally likely is the An = — 1 process. 

As an illustration, suppose we have e + e“ collisions at about a few hundred GeV, and we 
can tune the energy over a wide range. If the energy of the e + e“ pair does not cover a pass 
band, then the above ( B + L) violating processes will not happen. Let us assume that the e + e“ 
pair energy hits a pass band at energy E ~ 100 GeV (the second band in Table 1), which has 
a width of (using the more optimistic value in Table 1) 1CT 1 " TeV. Taking the energy spread 
of the electron beam to be 100 MeV, we see that only an exponentially small fraction (i.e., 
10— 173 ) of the e + e _ pair will lie inside the band. So such (B + L)-violating processes are indeed 
exponentially suppressed. 

The (B + L)-violating process in e + e“ collisions is very much suppressed due to the narrow 
band widths unless one can reach energies close to but below the sphaleron energy of 9 TeV. If 
one can reach that energy, say a couple of hundred GeVs below E sp h and has a narrow enough 
beam energy spread of about 1 GeV, one can tune the e + e“ energy to sweep over one band at 
a time to enhance the signal-to-background ratio. 
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Unfortunately, a 9 TeV e + e (or /i + /r ) collider is not feasible in the foreseeable future; so 
let us turn to proton-proton colliders, which offer a better chance of detection in the near future. 


7.2 Proton-Proton Collisions 

Here, we like to explore (B + L)-violating processes in proton-proton collisions based on the 
above Bloch wave analysis. In pp collisions, the fundamental processes are quark and/or gluon 
scatterings. As a quark’s momentum inside a proton has a wide parton distribution, the funda¬ 
mental parton-parton center of mass energy has a very wide spread relative to a typical band 
width. To enhance the cross-sections, let us consider the case when both partons are valence 
quarks. 

As an example, consider the following An = — 1 process, 

UlL + u 2 L -t e + /i + T + M 2 fc3CiC2S 3 (i3 + X, (7.2) 


where the subscripts are the color indices and X stands for all other particles. In this inclusive 
scattering, the spectator quarks from the two protons, namely u 2 d^ and u±ds, can combine with 
4 anti-quarks to form 4 mesons, leaving behind a single anti-baryon. 


For An = +1 process (1.4), consider, for example, 


U 1L + U 2L e /-t T b 1 b2t3CiC 2 C 3 UiU2U 3 diU2 + X, 


(7.3) 


a single pp collision can produce 5 baryons (plus baryon-anti-baryon pairs). Similarly, one can 
consider quark-gluon and gluon-gluon scatterings. Since gluons have no weak interactions, they 
must be converted to quarks or produce quarks for the (B + L )-violating processes to take place. 
Similar behavior can take place in an electron-proton collider. All these scatterings typically 
will produce additional massive bosons: Higgs Bosons, W ± and Z Bosons, photons as well as 
mesons. 


7.3 LHC Physics 

Let us consider LHC run at 8 TeV. The center-of-mass energy for any pair of quarks has energy 
E c < 8 TeV. For any of the above processes to be unsuppressed by an exponential tunneling 
suppression factor, it has to be inside a conducting band. For such a band with energy E less 
than 8 TeV, the band width A E <C 1CT 20 TeV. So at most a tiny fraction (of order 1CT 20 ) of the 
quark-quark collisions take place within such a band. Even when collision events cover many 
bands, we do not expect any observable signal, as the band width decreases exponentially fast 
as the band energy decreases. So not a single (B + L)-violating event is expected. 
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Now the situation can be very different in the LHC runs at 13 and 14 TeV. In the actual 
electroweak theory with sin 2 9 W = 0.23, the estimate is E sph = 9.0 TeV. Here, the gg scatterings 
can reach E(qq) > E sph = 9.0 TeV. A typical process can happen at energies covering a few dozen 
bands just below E sp h , where the band width varies from 10 GeV to MeVs. Therefore, there 
is a reasonable fraction of the scatterings taking place inside some bands, so (B + L)-violating 
events have a chance to be observed. Let us make a rough estimate of the event rate. 

Among parton-parton scatterings, only left-handed quarks have direct electroweak interac¬ 
tions. Let f(E(pp), E c (q L q L )) be the fraction of q L q L scattering with energy E(q L q L ) > E c (q L q L ) 
at proton-proton collision energy E(pp). Since the band width drops exponentially fast as 
E(QlQl) decreases, we may choose to introduce (rather arbitrarily) E^qr/jr,) ~ 9.0 TeV to sim¬ 
plify the estimates. Phenomenologically, with energies in units of TeVs, /(14,9) ~ 10 -6 to 10 -8 , 
based on the known parton distribution function for left-handed valence quarks |4l|443||. Among 
other parton-parton scatterings such as gluon-gluon, gluon-quark, g^g^j or q E qR scatterings, the 
partons have to convert to electroweak interacting particles first, which suppress their contri¬ 
bution to the (B + L)-violating process. The qq scatterings are expected to contribute little as 
well. Even in g L g L scatterings at energy above E c (q L q L ), radiation of a hard gluon can lower the 
energy of the qLqL P a i r to below E c (q L qL)- So we have to exclude the QCD g^g^, cross section 
as well. Let the fraction of g^g^ events that can participate in (B + L)-violating processes be 
Few ~ &ew(,QlQl)/&qcd((1l(1l)- We expect F EW to be relatively small. 

Even when a qLqL scattering can go through a ( B +L)-violating process without the tunneling 
suppression factor, it may choose not to. Suppose that cross section at energy E is a ew(E, An = 
0). Let the total g^g^ electroweak cross section at energy E be ctt(E) = cr EW (E, An = 0) + 
a(E,An = ±1), where we ignore the other (An ^ 0, ±1) contributions for simplification. Let 
the fraction of the (B + L)-violating processes among all electroweak processes be 

n(E) - a{E ' An = ±1) <- a{E ’ An = ±X) (7 4) 

°t(E) ~ o ew (E, An = 0) 


Now k = 0 within a band gap (no Bloch wave solution), while it is ’’tunneling” unsuppressed 
within a pass band, so we expect k < 1 at energies around and above the sphaleron energy. 


Averaging over a set of pass bands yields n < 1, even for energies above E sp h. To guesstimate 
the value of k, let us compare the (B + L)-violating process (1.3) to a An = 0 scattering 
with the same X as in process (1.3), i.e., q E q E —> qqX. We see that the extra leptons and 
quarks l e l p l r qqqqqqq in (1.3) have a threshold of about 20 GeV only, so phase space puts little 
constraint on the process (1.3) with energies above 9 TeV. Also, the (B + L)-violating part of the 
process (1.3) is non-perturbative; so, besides the tunneling suppression factor that disappears 
for E > E sp h, there are no powers of coupling factors that would have suppressed its rate. So we 
believe that k < 1, but not exponentially small. It will be very important to find out whether 
there are other factors that will make 1. 
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Let us now estimate the (B + L)-violating event number at LHC at E(pp) = 14 TeV. Note 
that there is a F 3 = 1/8 chance to observe 3 same sign charged leptons instead of one or more 
neutrinos in such processes. Let the total inelastic pp cross-section at 14 TeV be a(pp) ~ 80 
mb (millibarns) (see e.g., Ref |44j). Taking /(14TeV, 9TeV) ~ 10 8 (which is a crude but 
conservative estimate), and an integrated luminosity of L pp = 3000 /5” 1 (inverse femtobarns), 
we expect the number of (B + L)-violating events with 3 same sign charged leptons to be (units 
in barn and inverse barn) : 


Integrated event number 
~ cr(pp) ■ f (14,9) ■ F ew ■ k- F 3 - L pp 

~ (80 x 10~ 3 6)(10~ 8 )(10- 2 )(10“ 2 /3)(4)(3000 x lO 1 ^" 1 ) 

8 

~ 10 4 (7.5) 


where, for lack of better information, we naively take Few ~ 10~ 2 and the fraction n in Eq.(7.4) 
to be k ~ 10~ 2 /3. Crudely, we guestimate the number of such (B + L)-violating events to be 
10 4±2 in the LHC 14 TeV run. Clearly a better understanding of these and other possible factors 
is important. There is no question that increasing the pp energy E(pp) by just a few TeV above 
14 TeV will increase the event rate by more than an order of magnitude. 


Because the quark momenta are at the tail end of the parton distribution functions, we do 
expect /(14TeV, E c (qiqL)) to be substantially larger than /(13TeV, E c (qLqL))-, may be by as 
much as an order of magnitude. Other factors, such as the poorly understood F EW and n, 
probably vary much less as E(pp) increases from 13 TeV to 14 TeV; so they drop out in the 
ratio. As a result, the ratio of the rates of the (B + L )-violating processes may be close to 
/(14, 9)//(13, 9) ^ 1. Such processes may be detected by a careful comparison of the data from 
13 TeV with the data from 14 TeV. 


7.4 Remarks 


A couple of remarks on the phenomenology are in order : 


• As we go to higher energies in a proton-proton collider, larger An = — 1 — K processes 
become possible, 

qL + qL^r U p l T qqqqqqq + [UfJ T qqqqqqqqq ] A + V. 

For example, at E(pp ) = 100 TeV, the fraction f(E(pp ), \n\E sp h) for say |An| > 4 will have 
reasonable values for q E qi scatterings. In this case, a (B + L)-violating process that produces 
more than a dozen same sign leptons and a dozen of (anti-)b-quarks is quite feasible. 


One way to estimate the rate is via a sphaleron with CS number There are sphalerons 
with N = | where E(n/2) sp h > nE sp h. By replacing U(cp) (2.7) by U{rup ), the spherically 


27 






symmetric solution will have topological number | while its energy is greater than nE sp h . Lower 
energy multi-sphaleron solutions are in general only axially symmetric even without U(1)y 
{g' = 0) [45j : 46|. These multi-sphalerons are also repulsive at m# = 125 GeV, that is E(n/2) sp h > 
nE sp h and turning on g' does not change this feature. There is a direction in the field space 
where a periodic potential with these iV-sphalerons separating the vacua exists, so the Bloch 
wave physics should also apply. Typically, we expect such a process to become effective when 
the energy is close to or above that of a multi-sphaleron, E{n/2) sph > nE sph . To go from n = 0 
to n = 2 with some probability, the qpqp energy must be at least 18 TeV. 


• If a whole proton can fit inside a sphaleron, then we should consider the total proton- 
proton energy instead of the energy of its constituents in a {B + L)-violating process. However, 
the typical size of a sphaleron is m^ 1 while the size of a proton is Rp ~ 350m^, so sphaleron 
mediated baryon-number violating processes can take place only for the elementary constituents 
inside the proton, namely the left-handed quarks as described above. However, this point is not 
as obvious as it seems. 


As we have proposed that baryon-lepton number violating tunneling is not suppressed for 
some specific low energies, and as AKY shows that the size of the to-be-formed sphaleron (i.e., 
with 0 < n -C 1/2) at low energies can have a much bigger size, it is interesting to take a closer 
look on this issue. Following Eq.(5.6), we see that the size of gauge field part of a to-be-formed 
sphaleron is given by a -1 > , where a is given by Eq.(5.7). Demanding that a -1 > Rp , 

we find that n < 1CT 4 , which means the pp center-of-mass energy is below the lowest band 
energy, E < 4 GeV. So we cannot fit a proton into a to-be-formed sphaleron to generate a 
(B + L)-violating process. 


8 Summary and Remarks 


In this paper, we take advantage of the discrete symmetry of the periodic sphaleron potential 
to estimate the rate of the (B + L)-violating processes in the electroweak theory. The sphaleron 
energy E sp h = 9.11 TeV (or 9.0 TeV for phenomenology) measures the height of the potential 
barrier that separates vacua with different ( B + L) numbers. We write down the effective 
one-dimensional time-independent Schrodinger equation for this periodic potential where the 
coordinate is essentially g, which is closely related to the CS number n (see Eq.(3.1)). The 


choice of g yields a constant mass in the kinetic term. It is then straightforward to find the 
Bloch wave function and the conducting (pass) bands and their widths. When the quark-quark 
energy in a proton-proton collision is around or above E sp h, such ( B + L)-violating processes are 
no longer exponentially suppressed via tunneling. A crude estimate suggests that LHC at 14 
TeV has a good chance of detecting them, as some of these events will have 3 same sign leptons 
plus 3 b-quarks. 







Since the quark-quark energies close to but below E sp h are at the tail end of the parton 
distributions for proton-proton energies at 13 TeV, we expect a substantially higher rate of 
(.B + L)-violating processes at 14 TeV than at 13 TeV, so a comparison of the events at these 
two energies will be very interesting. If observed, they can provide a probe into other physics as 
well. For example, E sp h will be shifted if there are 2 Higgs doublets, or via quantum corrections 
or modified Higgs potential, topics not discussed in this paper. 


Bloch wave function is a special (particularly elegant) example of resonant tunneling: that 
is, tunneling can be very efficient when the resonance condition is satisfied. Resonant tunneling 
can happen only when there are 3 or more minima. This is very well understood in quantum 
mechanics, but much less so in quantum held theory. There is a very puzzling property in Helium- 
3 superfluidity, which has quite a number of distinct phases. If one starts in the A phase (while 
the actual ground state is in the B-phase) and tunes the temperature, pressure and/or magnetic 
held slightly, it will go to the B phase via a hrst order phase transition. As the held theory 
of He-3 superhuid is well studied, theoretical calculation shows that such a phase transition 
will take a minimum lO 20,000 years to happen 47 . Yet, in the laboratory, it happens readily. 
Now, being a superhuid, there are no impurities inside the sample; furthermore, the A—)-B phase 
transition must happen in the bulk of the sample away from the container wall, as wall effect 
stabilizes the He-3 to stay in the A phase. So there is nothing to seed nucleation bubble. Now, 
this puzzle may be explained by resonant tunneling 48 . Here, the (B + L)-violating process 


may be another case where we can directly check the resonant tunneling phenomenon in held 
theory in the laboratory. 


Our analysis, in particular the evaluation of the mass m, is based on two different approaches 
in the static approximation, namely that of the Manton and the AKY approaches. These two 
methods take different approximations to the actual problem. The differences in the band 
structures give us a rough idea of the uncertainties involved. Clearly a better evaluation is 
desirable. This probably requires a time-dependent treatment of the sphaleron potential and 
mass. 


Note that we have assumed the simplest Higgs potential in Eq. (2.1), with a quadratic plus a 
quartic term in <f>. Since only the Higgs vacuum expectation value and the Higgs mass have been 
measured so far, the Higgs potential may take other forms or there may be more than one Higgs 
doublet. This will change the sphaleron energy E sph (as well as the sphaleron potential shape) 
which in turn would change the threshold energy where the (B + L)-violating processes becomes 
unsuppressed. So the observation of such events at LHC will provide valuable information about 
the Higgs potential itself as well as the Higgs self-couplings. 


A better understanding of the (B + L)-violating processes at the colliders will surely improve 
our understanding of the (B + L)-violating processes during the electroweak phase transition in 
the early universe. We can also assess the possibility whether such processes will enhance the 
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anti-proton to proton ratio in our universe. 

A typical energy released in a nuclear reaction is of order MeV, while the energies released 
in the (B + L)-violating processes are of order of GeV (via the annihilation of baryons with the 
anti-baryons produced), roughly a factor of thousand bigger. Once we understand the physics 
better, we may find a way to take advantage of this release of energy. 
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A The function f(r) and h(r) 


We present in Fig. 9 the numerical results of the profile-function f(r) and h(r) of the gauge 
held and Higgs held respectively as defined in Eq.(2.7). 

Following the asymptotic expansion of /(r) and h(r ) at the origin, one can approximate 
them by 

f(r) ~ 1 — sech(amw^), h(r) ~ tanh (bmy/r) (A.l) 


and then determine a = 1.154 and b 
See Fig. 9 for comparison. 


1.056 by minimizing the energy Vm(I-i) (2.8) at ji — 7t/2. 
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Figure 9: (a) Numerical result (solid line) of f(r) and the approximate function (dashed line) 
f(r) « 1 — sech(1.154m^r). (b) Numerical result (solid line) of h(r) and the approximate 
function (dashed line) h(r) ~ tanh(1.056mvu^)- The numerical forms and the approximate 
functions almost overlap. 
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